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We introduce a new mechanism for generating magnetic fields in the recombination era. This 
Harrison-like mechanism utilizes vorticity in baryons that is sourced through the Bose-Einstein 
condensate of axions via gravitational interactions. The magnetic fields generated are on galactic 
scales ~ lOkpc and have a magnitude of the order of B ~ 10 -23 G today. The field has a greater 
magnitude than those generated from other mechanisms relying on second order perturbation theory, 
and is sufficient to provide a seed for battery mechanisms. 

PACS numbers: 95.35.+d, 98.80.-k 


Introduction — One of the most pressing problems in 
modern astronomy is determining the origin of magnetic 
fields in the Universe. Magnetic fields are observed on 
all scales, from small scales, such as inside our own solar 
system, to the largest bound structures, galaxy clusters 
[lJ-[sj]. In fact, recent observations have even detected 
an intergalactic magnetic field existing in the void re¬ 
gions of the cosmic large scale structure, with magnitude 
Higmf ~ 10 -18 — 10 _15 G (e.g., Refs. 0,0]). Despite 
their prevalence, there is a large amount of uncertainty 
in how the first magnetic fields were created. 

Setting aside the possibility that magnetic fields were 
present as initial conditions which is incredibly unappeal¬ 
ing, we are yet to fully understand the origin of the first 
fields. These primordial seeds need only be small as there 
are several ways in which they can then amplified by as- 
trophysical processes; for example, adiabatic contraction 
or turbulent dynamos during structure formation (see, 
e.g. Refs. [6)-|8j] for a review). The seeds are required 
to have an amplitude in the range of 10~ 3 ° — 10~ 2 ° G, 
with the specific magnitude depending on the details of 
the dynamo model. There has been much work over the 
years addressing the generation of the primordial seed 
field. 

Harrison 1,0 was one of the first to attempt to ex¬ 
plain the origin of the seed magnetic field generated by 
vorticity in a rotating protogalaxy prior to decoupling. 
This provided a seed field of the order of B ~ 10 -19 G, 
which is large enough to source galactic dynamo mech¬ 
anisms which enhance this initial seed field to currently 
observed magnitudes. However, since (at linear order in 
perturbation theory) vorticity decays [llj, there is no way 
to support the vorticity in a post-recombination universe, 
and so this mechanism of magnetic field generation was 
criticized [12]. Later, Mishustin & Ruzmaikin [13| in¬ 


vestigated the generation of magnetic fields in the post¬ 
recombination plasma, finding fields with magnitude of 
around 10 _17 G (evaluated at z ~ 100). Similarly to 
Harrison, this work required the existence of primordial 
vorticity. Another recent piece of work used vorticity 
from the texture scenario of large scale structure forma¬ 


tion [14]. However, the resulting field is too weak to act 
as the primordial seed. The other important seed field 
mechanism in the early universe is the Biermann battery, 
either due to shocks 151 or during reionization EE- 

Other mechanisms for generating magnetic fields in the 
very early universe have been studied comprehensively 
in the literature. These roughly include fields generated 
during inflation, with a breaking of the conformal invari¬ 
ance of electromagnetism jl8| . during phase-transitions 
00 or during (p)reheating 0,|22|. These all have 
issues, and sustaining magnetic fields in the very early 
universe proves to be difficult. 

One interesting method for generating magnetic fields 
around recombination builds upon the ideas of Harrison, 
and uses second order cosmological perturbation theory 
00. While vorticity decays at linear order in per¬ 


turbation theory, there are source terms at second order 
- these look very much like the baroclinic term in the 
Biermann battery - that allow for vorticity generation 
26- 281. It is therefore possible that this vorticity comes 


hand-in-hand with a magnetic field. There have been nu¬ 
merous works to this end 0-0], that all obtain fields 
with roughly the same magnitude of 10~ 26 — 1CU 23 G at 
recombination. 

In this paper, we consider a new method in which 
the vorticity can grow around decoupling by consider¬ 
ing axion dark matter. The axion is a viable dark matter 
candidate [35 37J, alongside the other major candidate 
- the broad class of weakly interacting massive particles 
(WIMPs). Although axions and WIMPs are similar, the 
axion is truly a quantum scalar field 0], and therefore 
it is expected that there will be observational differences 
between the two. In particular, it was recently realized 
that axion dark matter can form a Bose-Einstein con¬ 
densate (BEC) 0 , 0 . Furthermore, the axion BEC 
can exhibit vortices, at the zeroes of the wave function 


41]. It is this vorticity that we will use in the present 


work to generate magnetic fields in the early universe. 

The axion Bose-Einstein condensate — The ax¬ 
ion is motivated by the Peccei-Quinn solution of the 
strong CP problem, and its mass is thought to be around 
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10 -5 eV/c 2 . Cold axions are one of the leading dark mat¬ 
ter candidates. Recently it was shown that axions form a 
re-thermalizing Bose-Einstein condensate (BEC) through 
gravitational self-interactions when the photon temper¬ 
ature was around 500 eV [H, [Z(l 0 ■ The axion BEC 
interacts gravitationally with baryons with a relaxation 
rate 


Eq ~ 47t Gnmrrib —— 
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where n and m are the number density and mass of ax¬ 
ions respectively, m b is the mass of baryons which is of 
order 1 GeV, £ ~ is the correlation length of the axion 
BEC, A pb ~ i/3m b T is the momentum dispersion of the 
baryons, where T is the photon temperature. 

For T < 1 keV the dominant interaction between pho¬ 
tons and baryons is Compton scattering. The relaxation 
rate for an electron to gain or lose energy by Compton 
scatte ring off photons is known from standard cosmology 
to be 431 
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where T 7o is the present day photon temperature and 
f \hr is the present day physical baryon density param¬ 
eter. It can be shown that both the ratios Tq/H and 
r g/T e are greater than one around matter radiation 
equality and keeps increasing thereafter. Therefore the 
baryons thermalize with the axion BEC. 

Vorticity from tidal torque — In the standard pic¬ 
ture of structure formation baryons collapse onto dark 
matter overdensities. Tidal torque from nearby inhomo¬ 
geneities imparts angular momentum onto such proto¬ 
galaxies. It was shown in Ref. 41] that when a system 
of axion BEC acquires angular momentum, the axions 
thermalize and most of them go to a state with mini¬ 
mum rji = ei — £(t)li, where e,; and U are the energy and 
angular momentum respectively of the i th state and £(i) 
is the angular velocity of the system which grows with 
time as angular momentum grows. This lowest y state 
has non-zero vorticity and therefore the axions acquire a 
net rotational velocity field. 

The baryons being in thermal contact with the axion 
BEC are dragged along with the axion flow and acquire 
the same rotational velocity field as the axions. The 
baryons therefore acquire vorticity from tidal torquing 
as a result of thermalization with the axion BEC. It 
should be noted that in general tidal torque on collision¬ 
less particles cannot generate rotational flow 44]. Before 
shell crossing the baryons behave like collisionless par¬ 
ticles since dissipative processes and shocks are absent. 
Therefore if dark matter is made of only WIMPs then 
there is no vorticity. 

Magnetic fields — We will now show how vorticity 
in the baryon fluid can generate a magnetic field. We 


follow an approach similar to that of Ref. [13j. As re¬ 
combination begins at redshift z T ~ 1500, the photon 
free streaming length grows rapidly such that, on galac¬ 
tic scales, it can be treated as a homogeneous radiation 
background. The charged particles moving in this radi¬ 
ation background experience a Thomson drag force, 451 
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v, where <tt is the Thomson cross-section, p 7 


is the energy density of photons and v is the velocity of 
the charged particle relative to the background radiation. 
Because electrons are much lighter than protons, their 
acceleration due to this force is much greater than that 
of protons. Neglecting the effects from neutral species, 
the equations of motion for electrons and protons with 
velocities v e and v p respectively are 
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where r ep is the characteristic time for momentum ex¬ 
change via Coulomb scattering, m e and m p are the 
masses of the electron and proton, respectively, e is the 
charge of the electron, and a grav is the acceleration due 
to gravitational interactions with the axion BEC. Since 
we are interested in showing how the vorticity in baryons 
sourced by the axion BEC can generate magnetic fields, 
we have neglected the electron and proton pressure terms 
in the above equations. 

The current density is defined as J = en e {y p — v e ), 
where n e ~ n p , by local charge neutrality. The Thomson 
drag term becomes negligibly small after z ~ 900 when 
the timescale of Thomson scattering is greater than the 
Hubble time. We have assumed that initially there is no 
magnetic field and neglected the back reaction from the 
generated magnetic field. 

On taking the difference of Eqs. © and © and ne¬ 
glecting the Thomson drag term for protons, followed by 
the curl, we arrive at an equation for the vorticity of 
electron fluid, u5 e , 




—V x E 

m e 


, 4<t t p 1 _ 
+ w 
6m e c 


e 


— 2 V x 


Je 

m e cr 


( 5 ) 

where a is the conductivity of the background medium. 
The LHS of Eq. © is negligibly small compared to the 
first term on the RHS on galactic scales [l|. The last 
term on the RHS of Eq. © is proportional to the mag¬ 
netic diffusion term which can be neglected because of 
the high conductivity of the background medium. We 
are therefore left with 
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On invoking the Maxwell equation 


1 dB 
c dt 


-VxE, 


Eq. ©, becomes 
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Of course, these calculations are performed in a static 
universe, therefore we must transform to the expanding 
universe in which we live. On doing so, Eq. © becomes 


1 d(a 2 B) _ Aa T p l0 a 4 _ 
a 2 dt ~ 3e We( ' 
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where a(t) is the scale factor and a subscript zero denotes 
the present-day value of a quantity. 

Let us consider a galaxy sized (~ 10 kpc) spherical 
overdensity of axion BEC onto which baryons are falling. 
Tidal torque imparts the same specific angular momen¬ 
tum to the infalling matter. Thermalization with the ax¬ 
ion BEC results in vorticity in the baryons which is of the 
order ui ~ L/MR 2 , where L is the total angular momen¬ 
tum, M is the total mass of the infalling baryons and R 
is the size of the protogalaxy. Following Peebles 46[ , the 
angular momentum of a protogalaxy grows as t 5 / 3 in the 
linear regime, which implies that the vorticity grows as 
f 1 / 3 . At z ~ 10 the protogalaxies reach their turnaround 
radius after which they begin to collapse and separate 
from the background. We denote this redshift by z* in 
the following. From this time onwards the evolution is 
complicated to handle analytically as non linear effects 
play a significant role. We make an estimate by con¬ 
sidering that the angular momentum of the protogalaxy 
is conserved per comoving volume after they separated 
from the background, so the vorticity decays like t -4 / 3 . 

To summarize in terms of redshift we have 


f w 0 (1 + z*) 5/2 (l + z) 1/2 , z* < z < z r 
uj[z) = < 

|wo(l + z) 2 , 0 < z < z* , 

( 10 ) 

where u >o is the present day value of the vorticity which, 
for our galaxy, is ujq ~ 10“ 15 s _1 . Expressing Eq. © 
in terms of redshift and using the above expression for 
vorticity we get an equation which can be integrated from 
the beginning of recombination upto z ~ 900 when the 
battery shuts down. For z r > z > 900, we have 
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The magnetic field grows up to z ~ 900 when it 
has magnitude B ~ 10 -1 ' G. After this time it is 
frozen into the residual free charges and decays with the 
expansion of the universe. The magnetic field today has 
a magnitude of Bq ~ 10 -23 G on scales of order 10 kpc. 


Discussion — In this paper, we have investigated the 
generation of magnetic fields from vorticity in the recom¬ 
bination era. We have used a Harrison-like mechanism, 
with the novelty lying in the fact that the vorticity is 
not assumed, but rather is inherent in the Bose-Einstein 
condensate of axions. This provides a natural source of 
vorticity which is present only for axion dark matter. The 
magnetic field sourced by this vorticity has a magnitude 
of B ~ 10 —1 ' G peaking at redshift z = 900, on scales 
of 10 kpc whose value today is of order 10 -23 G. The 
magnetic field generated through this process acts as a 
seed for astrophysical amplification mechanisms through 
the later stages of galaxy formation. There are several 
different dynamic mechanisms which can amplify seeds 
by upwards of ten orders of magnitude 0E3 , and result 
in the observed fields of the order of a few microGauss 
at redshift less than one. 

Furthermore, the magnetic field generated from axion 
dark matter is larger in magnitude that those created by 
mechanisms relying on higher order fluctuations within 
the standard ACDM cosmological model. Therefore, this 
allows for less effective amplification mechanisms to en¬ 
hance the primordial seed to the observable size. 

Finally, we should note that taking into account 
effects on how the baryons collapse more than the 
dark matter halo (e.g. Ref. [48j]), the field could 
be diluted by a factor of (20) 2 in the inter galac¬ 
tic medium (IGM) compared to the disk. This will 
result in a field with magnitude B ~ 10 _25 G in the IGM. 
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